Supercurrent in Long SFFS Junctions with Antiparallel Domain Configuration 
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We calculate the current-phase relation of a long Josephson junction consisting of two ferromag- 
netic domains with equal, but opposite magnetization h, sandwiched between two superconductors. 
In the clean limit, the current-phase relation is obtained with the help of Eilenberger equation. In 
general, the supercurrent oscillations are non-sinusoidal and their amplitude decays algebraically 
when the exchange field is increased. If the two domains have the same size, the amplitude is in- 
dependent of h, due to an exact cancellation of the phases acquired in each ferromagnetic domain. 
These results change drastically in the presence of disorder. We explicitly study two cases: Fluc- 
tuations of the domain size (in the framework of the Eilenberger equation) and impurity scattering 
(using the Usadel equation). In both cases, the current-phase relation becomes sinusoidal and the 
amplitude of the supercurrent oscillations is exponentially suppressed with h, even if the domains 
are identical on average. 

PACS numbers: 74.45.+c,74.50.+r 



I. INTRODUCTION 

Hybrid systems containing superconducting and ferro- 
magnetic elements gained recently a lot of attention due 
to experimental progress as well as possible applications 
in magnetoelectronics and quantum information. The- 
oretical studies are revealing a variety of new features, 
making these system generators of novel theoretical con- 
cepts. 

It is a common knowledge that current in hybrid nor- 
mal metal - superconductor (NS) systems flows by means 
of Andreev reflections: an electron in N is reflected from 
the NS interface as a hole with the opposite charge and 
velocity, fmagine first that the piece of normal metal is 
ballistic. An electron at the Fermi surface is reflected 
as a hole at the Fermi surface, and they propagate in 
the normal metal with the same phase. If the electron 
is taken at a finite energy E (counted from the Fermi 
surface), a momentum mismatch 5p = 2E/vp between 
this electron and the reflected hole appears, vp being the 
Fermi velocity. 

Consider now an interface between an (s-wave) super- 
conductor and a ferromagnet. Electron and hole have 
opposite spin directions, and the exchange field h in the 
ferromagnet leads to a Zeeman splitting of energies of 
the two different spin projections. Thus, even an elec- 
tron and a hole at the Fermi surface acquire the momen- 
tum mismatch 2h/vF] hence their relative phase grows as 
8ip = 2hx/(hvp), where x is the distance from the inter- 
face. This affects phase-sensitive physical quantities like 
the supercurrent in superconductor - ferromagnet - su- 
perconductor (SFS) junctions: It becomes an oscillating 
function of the thickness d of the ferromagnetic layer, 
with a period hvp/2h. If, furthermore, the ferromag- 
net is diffusive, the oscillating behavior is accompanied 
by an exponential decay cx exp(— (h/hD) 1 /' 2 ^, where D 
is the diffusion coefficient. Typically, h is much larger 



than the superconducting gap A, and thus the length 
scales related to the magnetic field are much shorter 
than the superconducting coherence length £, Hvf/ A and 
(HD/A) 1 ' 2 in the clean and diffusive case, respectively. 
In other words, the proximity effect is suppressed in the 
ferromagnet. 

This qualitative discussion suggests that the main ef- 
fect observed in SFS contacts is oscillations of the su- 
percurrent with the thickness of the ferromagnetic layer 
— the transition from a so-called 0-state (energy of the 
contact is minimum for zero phase difference between the 
superconductors) to a 7r-state (energy is minimum for a 
phase difference it) . This topic was at the focus of atten- 
tion since the early exploration of the field [l| . Theoreti- 
cally, the 7r-state was predicted in a variety of SFS junc- 
tions: Ballistic 0,11110, short diffusive [10, long dif- 
fusive Pa, IE a, ferromagnetic insulating barrier ffl Q IIlCl . 
ballistic U^M, and diffusive El H Eflllllll 
junctions with a barrier separating two ferromagnetic lay- 
ers, and ballistic \Vt^ and diffusive 0] with two tunnel 
barriers. The transition to the 7r-state was recently ob- 
served experimentally in SFS junctions EHl El \2 ^ . 
All these observations are limited to small thickness of 
ferromagnetic layer (s), d < (hD/h) 1 / 2 . For thicker lay- 
ers, supercurrent does not exist. 

In this situation, it is useful to understand how one 
can enhance the proximity effect. Several options have 
been recently discussed in the literature. First, the above 
qualitative argument assumes that the pairing between 
an electron and a hole participating in the Andreev re- 
flection is singlet - they have opposite spin projections. 
Obviously, if the superconductor allows for a non-trivial 
symmetry of the order parameter, this needs not be the 
case, and triplet pairing between an electron and a hole 
with the same spin projection can arise. Since triplet- 
paired electron and hole at the Fermi surface have no 
momentum difference, they can propagate with the same 
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phase and enhance the proximity effect. Coupling of 
two d-wave superconductors via a ferromagnetic layer has 
been considered in Ref.|3 Moreover, triplet pairing can 
even appear in a contact of an s-wave superconductor and 
a ferromagnet, pro vided the magnetization in the latter 
is non-uniform |25l l2fil l27l | . In this case, the proximity 
effect survives at the same distance £ from the interface 
as in non-magnetic metals. Indeed, the supercurrent in 
SFS junctions with non-uniform magnetization is con- 
siderably enhanced 0- We also mention that the super- 
current in a long diffusive SFS junction is exponentially 
suppressed only on average; phenomena related to the 
proximity effect still occur in such a junction as a result 
of mesoscopic fluctuations around average quantities (2^| . 
Finally, if the ferromagnetic layer is split into domains, 
the coherence can be preserved if an electron and a hole 
propagate between the superconducting electrodes along 
the two sides of a domain wall ^3 ■ 



In this Article, we explore a different way to enhance 
the supercurrent in SFS junctions. Imagine first that the 
junction is ballistic and the ferromagnetic layer consists 
of two domains with opposite directions of the magnetiza- 
tion, as shown in Fig. 1. Triplet pairing is not generated 
in this geometry. Consider an electron and an Andreev- 
reflected hole propagating from left to right between the 
superconducting electrodes. They first acquire the rel- 
ative phase Sifi = 2hx\/ '(Hvf), x\ being the distance 
traversed in the first ferromagnetic layer. However, in 
the second layer the exchange field has the opposite sign, 
and the phase gain 6(f2 — —Ihx^j (Jvvf) partially com- 
pensates 5ip±. For x\ — X2 we have full compensation: 
The ferromagnetic bilayer behaves as a piece of normal 
(not ferromagnetic) metal, and the proximity effect is 
fully restored. Indeed, previous studies of SFS contacts 
where the two ferromagnetic domains were separated by 
a barrier, found that the supercurrent in the antiparal- 
lcl domain confi gura tion is enhanced with respect to the 
parallel one [Til 1131 ITU ITsj . If the domains are identi- 
cal, there is no transition to the 7r-sate in the antiparallcl 
configuration. 



Below, we consider such a situation quantitatively. 
Section [H] treats a ballistic SFFS junction with two fer- 
romagnetic domains parallel to the superconducting in- 
terfaces. We show that this system behaves as a bal- 
listic SFS junction with an effective exchange field. If 
the widths of the two domains are the same, this effec- 
tive field vanishes. In the next two Sections, we study 
the effect of disorder in the same system and show that 
supercurrent in diffusive SFFS junctions decays exponen- 
tially with their width, similarly to SFS contacts without 
domains. We consider long junctions, d 3> £, and assume 
that the superconducting electrodes do not influence the 
magnetic structure of the contact. 



II. CLEAN SFFS CONTACT 

We consider first a system of two clean ferromagnetic 
strips |30| with the thicknesses d\ and di and antipar- 
allel orientations located between two superconductors 
(Fig. 1). The dynamics of quasiparticles in this system 
are described by the Eilenbcrger equation, 



WFTiVg^ir, n) = (iuj T ha)f 3 + A, g a (r, n) 



(1) 



Here the semi-classical Green's function g a is a matrix in 
Nambu space, 



9o 



9a fa 
fa ~9a 



which describes the singlet pairing (the triplet component 
is not generated in our geometry), and the spin index 
fj = ±1. The exchange field h is zero in the supercon- 
ducting banks, and has antiparallel orientations in the 
ferromagnets: The upper/lower signs in Eq. corre- 
sponds to the left/right ferromagnet [h > 0). To stay 
in the framework of the semi-classical consideration, we 
have assumed that the Zeeman splitting h is much weaker 
than the Fermi energy, but can be arbitrary in compar- 
ison with the superconducting gap A. We put h = 1; it 
will be restored in the final results. 




FIG. 1: Setup with two ferromagnetic domains with antipar- 
allel configurations. 



In this Article, we consider the case of a long contact: 
The thicknesses of both ferromagnetic layers are much 
larger than the superconducting coherence length, c?i.2 2> 
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Hvf/A. Then the matrix A can be taken in a piecewise 
approximation: It is zero in both ferromagnets, and 



Ae tx 
-Ae- lx 



in the superconductors. Here \ = —Lp/2 and x — y/2 in 
the left and the right superconducting bank, respectively. 

In the bulk superconductor far from the contacts the 
Green's function is isotropic and equals for |o>| < A 



bulk 



VA 5 



LO 



-iAe lx 

— LO 



(2) 



In addition, the Green's function and its derivative must 
be continuous at each interface. 

We introduce the coordinate x parallel to n and di- 
rected from left to right. Let us choose x — at the 
boundary of the left superconductor; then x = d\j cos 
at the interface of the two ferromagnets, and x = (d\ + 
d 2 )/ cos 9 at the boundary of the right superconductor. 
The quasiparticles in the clean system move along a 
straight line (Fig. 0. It follows from Eq. that the 
normal component g a (r,n) is constant along the trajec- 
tory inside the ferromagnets. The calculation gives 



, . VA 2 + lo 2 sin a + + cos a) 

5a(n) = — 1 ■ /A2^ fTT T s> ^ 

lo sin a + i\J lV- + lo z (1 + cos a) 

where the phase a accumulated along the trajectory is 

2ilo di + d 2 2ha di — d 2 

a = ! n x tp, n x = ±l. (4) 

vp cost) vp cos 6 

The supercurrent density is expressed as follows, 



J 



-nrevp 



V ^Z T ^2 J dn 9nn, 



(5) 



where v is the density of states. For h = Eq. (JjjJ gives 
the supercurrent of a long clean SNS (non-ferromagnetic) 
junction, as considered in Ref. l3ll which we follow in 
the general case. The expression is even in lo; for 
zero temperature (case of interest here) the summa- 
tion can be replaced by an integration over frequencies. 
We subsequently introduce a new integration variable 
lo = A sinh u and arrive at the intermediate expression 



/>OG /* 7r /2 

j = 2et>i?£ / A^^ / ducoshu / d9 cos 
„ JO Jo 



(6) 



x Im tanh 



. A . , d x +d 2 d x -d 2 op 

u + A sinh u + ilia + i— 

vfcosO vfcosO 2 



For long contacts, Adi y2 ^S> hvF, the first term in the 
argument of the hyperbolic tangent can be disregarded. 
Using the identity 



Imtanhj/ = 2^(— l) fe Im e 



2ky 



we obtain the final expression for the supercurrent, 

■ sin kip 



J 



4ev 2 F vh^ (-l) fe+1 



dl + ^ 2 k~i k 



dx 2kh(d\ — d 2 )x 

cos ; . (7) 



/i x 3 Vx 2 - 1 v F h 
For h = 0, we return to the clean long SNS contact 



3 = Jo 



(-1) 



fe+i 



• sin kip, 



(8) 



fe=i 



where jo = irev F vh/(d\ + d 2 ). This describes the well- 
known sawtooth current-phase relation found earlier in 
Ref.|H 

For strong magnetic fields, h ^> hv F /\di — d 2 \, the in- 
tegral over dx in Eq. 0, which corresponds to summing 
over all possible trajectories in the ferromagnets can be 
calculated in the saddle-point approximation. As a result 
we find the current-phase relation 



2.7'cm 



X cos 



Vf^ 



ith\di - 
'2kh\di 



fe+i 



vfIt- 



1 fc=i 

4 



fe3/2 

sin kip. 



(9) 



k=l 



We note, first of all, that the amplitude of the supercur- 
rent oscillations as a function of ip decreases algebraically 
with the exchange field, as \Jh,VF jh\d\ — d 2 \. This is a 
direct consequence of the fact that we summed over all 
possible trajectories, and hence averaged over the differ- 
ent phases acquired during propagation in the ferromag- 
netic domains along these trajectories. Secondly, as far 
as the phase- dependence of the supercurrent is concerned, 
it is in general neither sinusoidal, nor saw-tooth-like. In 
Fig. [21 we plot j(tp) for various values of h\d\— d 2 \/hvF ~ 
10, such that the saddle-point approximation is reason- 
able. We see that, as a function of the exchange field, 
the supercurrent changes sign at a given phase differ- 
ence. Thus, depending on the parameter h\di — d 2 \/hvp, 
the junction either favors a 0-state or a 7r-state. We fi- 
nally note that for d 2 = 0, Eqs. J7J and © give the 
supercurrent for a (single-domain) clean long SFS junc- 
tion. This is, to our knowledge, a new result as well. It 
implies in particular that a clean SFS junction can also 
be a 7r-junction, in accordance with previous results for 
different types of SFS hybrid structures. 

The important conclusion for the general case is that 
for a two-domain contact the result is exactly the same 
as for a SFS junction with the thickness d\ + d 2 and 
the effective exchange field h e f = h\di — d 2 \/(d\ + d 2 ). 
In particular, if the thicknesses are the same, d\ = d 2 , 
the magnetic field drops out - we obtain the sawtooth 
current-phase relation (JSJ) like for a SNS contact. In the 
language of Eilenberger equations, this statement is obvi- 
ous: Indeed, the only quantity sensitive to the magnetic 
field is the phase a accumulated along the trajectory. 
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FIG. 2: Supercurrent j in units of jo as a function of ip for 
various values of h: h\di — di\jhMF = 10.0 (a), 12.0 (b), and 
14.0 (c). 



Since each trajectory is a straight line, each layer con- 
tributes with a weight proportional to its thickness and 
with the sign depending on the direction of the exchange 
field. This result is readily generalized to the case of 
many ferromagnetic layers in the antiparallel configura- 
tion 1331. 



III. DISORDER AVERAGING 

Now we discuss how our two main observations for the 
supercurrent — power-law decay with magnetic field and 
the independence on the magnetic field in the symmetric 
case d\ = d 2 — react to the presence of disorder. Before 
performing this difficult task in Sect ion ITVl bv solving the 
Usadel equations, we try to use an easy way to under- 
stand the effect of impurities in this Section. We intro- 
duce randomness in the thicknesses of the layers (surface 
randomness). This simple and transparent calculation 
provides us with results which are clear qualitative pre- 
dictions to be compared with the conclusions extracted 
from a more complicated analysis of the Usadel equa- 
tions. 

We start from Eq. (JJJ , and imagine that the interfaces 
as presented in Fig. ^ are not straight, but exhibit small 
fluctuations in position. Since there is no scattering at 
the interfaces, the only effect of such fluctuations is that 
the thicknesses of the layers become random variables, 
and the supercurrent Q must be averaged with respect 
to this randomness. Let us take a Gaussian distribution 
for the difference d\ — d 2 , 



P(di 



exp 



di + d 2 f 



(10) 



where a <C di,d2 has the meaning of a typical scale of 
the interface fluctuations, and di are the averaged values 



both thicknesses. Averaging Eq. J7J), we obtain 

) k+i . , r dx 

■ y - — ; sin kip j 

7T 



4j ^ (-!)» . , 

k=l 



x exp 



k 

k 2 h 2 a 2 x 2 

v%h 2 



c 3 ya 



1 



cos ■ 



2kh(d\ — d 2 )x 



(11) 



where jo = 7rev^ph/_(di + d 2 ). In strong fields, h 3> 
hvp/a 3> hvp/(di — d 2 ), the integral is calculated in the 
saddle-point approximation, and only the term with k = 
1 survives, 



y nh\di — d 2 
(2h{d 1 - d 2 ) 



vph 



exp 

7T 

" 4 



h 2 a 2 
vfE 2 



sm if. 



(12) 



Thus, the averaging procedure brings out two, qualita- 
tively new features: (i) at high fields, the current-phase 
relation becomes sinusoidal; (ii) the amplitude of the su- 
percurrent oscillations decays exponentially, rather than 
algebraically with h. In addition, the exchange field still 
modulates the phase of the oscillations, and can drive the 
contact to a 7r-state. The property (i) stems from phase- 
averaging over diffusive trajectories and is a common fea- 
ture of all long disordered SNS junctions (cf Ref . I34T) . 

Eq. (|12|) does not apply to the symmetric case d\ — d 2 . 
In this situation, for h vf/cl, we have 



J 



2 jo VF^ 
\/7T ha 



exp 



sm ip. 



(13) 



We see that even for the symmetric case the exponen- 
tial dependence on magnetic field persists. It reflects the 
fact that a quasiparticle moving along a single trajec- 
tory spends in general unequal times in both layers, and 
thus the contribution of each trajectory is magnetic field 
dependent. However, there is no additional oscillating 
factor due to the magnetic field: a symmetric junction 
is never in the 7r-state. These features are confirmed 
qualitatively in the next Section, where we analyze the 
behavior of a symmetric diffusive SFFS junction, using 
the Usadel equations. 



IV. DISORDERED SFFS CONTACT FROM 
USADEL EQUATIONS 

We now consider a diffusive SFFS junction in the sym- 
metric case d\ = d 2 = d/2. The junction is again as- 
sumed to be long, d ^> (hD/ A) 1 / 2 , with D being the 
diffusion coefficient. 

If the exchange magnetic energy does not exceed the in- 
verse elastic scattering time, h <C h/r, the Green's func- 
tion is almost isotropic, and the system can be described 
by Usadel equations, 

Dd z [F a d z F+ - F+d z F a ] - 2(AF+ - A*F a ), (14) 

Dd z [G a d z F a - F a d z G a ] = 2{lo ± iha)F a - 2AG a , 
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with the constraint G 2 + F a F+ = 1. Here, as usual [34| . 



and it actually only depends on the distance z from the 
ferromagnet-ferromagnet interface (Fig. The up- 

per/lower signs describe the regions —d/2 < z < and 
< z < d/2, respectively, and A ~ iAexp(i\) in the 
superconductors. In the following, we suppress the spin 
index a where it does not lead to ambiguities. 

Following Ref. 1.3 ll . we solve the constraint by introduc- 
ing two complex- valued fields 9 and 77, 

G = cos0, F = sm9e lv , F+ = sintfe""'. (15) 

The equation for rj in the ferromagnets becomes 

5 2 (r/sin 2 6») = 0, (16) 

with the boundary conditions r)(±d/2) = z f(p/2. The 
first integral yields 



V 



I 



(17) 



where / is an unknown constant. The current is ex- 
pressed via this constant, 



(18) 



To ensure the current conservation, / must be the same 
in both ferromagnetic layers. It is important, however, 
that we do not assume that the current is conserving - it 
follows naturally from the consistency of our solution. 

Using Eq. (|16|) . we also write the equation for 9 in 
ferromagnets, 



D9" = DI 



COS! 



sin 



2(lj ± iha) sin#, 



(19) 



with the first integral 



DO' 



DI 2 



A(uj ± iha) sin 9 + const. (20) 



Now, for the long junctions, the boundary conditions for 
9 at z = ±d/2 are essentially the same they would be 
at the interface between a semi-infinite superconductor 
and a semi- infinite ferromagnet. To find these boundary 
conditions, we write the corresponding equation for the 
superconductors, 

D9' 2 = -4ui sin 9 - 4A cos 9 + const, 

Taking into account that in the bulk superconductor 9 = 
7r/2, in the bulk ferromagnet 9 = 0, and requiring the 



continuity of 9 and 8' at the interface, we obtain the 
following boundary conditions, 



u ± iha{l - cos 6) + A(l - sin ( 



sin 2 9 



= 0, 



at z = ^d/2. Although our equations describe the behav- 
ior of an SFFS junction for an arbitrary relation between 
h and A, we concentrate in the following on the case 
T, h <g; A. As we show below, in this situation the cur- 
rent / is exponentially small, and the boundary condition 
for 9 reduce to 9(z = ±d/2) = tt/2. 

Since the Usadel equations posses obvious symmetries 
9 a (Lo) = 9- a (uj) + 7r, rj a {uj) = r/-o-(oj) + n, in the sequel 
we only consider lo > 0. 

The field 9 must rapidly decay away from supercon- 
ductors and stay exponentially small within the ferro- 
magnets. We start first solving Eq. I|2U|I at z <§; d, where 
9 <C 1, and the trigonometric functions can be expanded. 
Then Eq. I|2U|I can be integrated. The solution is too 
cumbersome to be written down here, its asymptotics for 
\z\ — > 00 are 




with the notations 9$ = 9(z = 0), 7 = 8'(z = 0), and 



a = -7= \/ V 'oj 2 + h 2 + lu; f3 = ±-7=V V 'u 2 + h? — uj. 
v 2 v 2 

Next, we solve Eq. (I20|l close to the interfaces, \z — 
d/2\ <C d. We assume that I/80, 7 are both exponentially 
small (to be checked later) and obtain 

9 it I [2 , „ ... , d.\ 

tan - = tan-exp l^j —(a + if3a){\z\ - -) I . (22) 

Far from the interface, 8 <C 1, the solution becomes 
exponential. Matching the exponential asymptotics of 
Eqs. (EU and (E3), we find the condition 




D7 2 



2{u> ± iha) 



DF 2 



64 tan — exp 



29 2 {u±iha) 
-(a + i/3a)d I . 



(23) 



We now integrate Eq. (|17fl . Since 9(x) grows exponen- 
tially away from x = 0, the sine in the denominator can 
be replaced by its argument. We then find 



DP 



2(u) ± iha) 



D7 2 



2(u J + l ha)) tan \2 ±m )' 
(24) 
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with i]o = 77(0) . We proceed by calculating the four quan- 
tities /, 9q, an d 7- The result is 




(25) 



Note that / docs not depend on a. It can be easily 
checked that I/9o and 7 are exponentially small, which 
justifies the approximations we have made to arrive at 
Eq. 

Now we calculate the supercurrent according to 
Eq. (fT5|) . For high temperatures T ^> HD/d 2 ,h only 
the term with lo = 7rT is important, and we obtain 



= v^io.diff 



irkBTcP 

hD 



3/2 



x exp \ --^=\[txT + yjh 2 + tt 2 T 2 ^ siny, (26) 

where we introduced jo,diff = 128(\/2— l) 2 evHD 2 /d 3 . In 
high magnetic fields h^> T, hD /d 2 the terms with uo < h 
contribute, 



J — Jo.diff 



hd 2 \ 3 ' 2 ( f~h\ . 

— — exp —\ a sin 05. 



(27) 



We note the two main features of the solution in the diffu- 
sive case. First, the current-phase relation is sinusoidal. 
This corresponds to the result for the long diffusive SNS 
contact j34|. Then, the supercurrent decays exponen- 
tially with magnetic field, in contrast to the power-law 
decay in the clean case. 

Similarly, we can treat a single-layer SFS junction of a 
thickness d. The result for h HD/d 2 ,T reads 



V. DISCUSSION 

We considered the behavior of the supercurrent in 
long SFS junctions. We obtained new expressions for 
single-domain ballistic and diffusive contacts and con- 
firmed that the to 7r transition can be induced in these 
systems. However, our main focus is on the situation 
when the ferromagnetic region is split into two ferromag- 
netic domains with equal but opposite magnetization. 
In the ballistic case, this system behaves as a single- 
domain SFS junction, with the effective exchange field 
h e t — h\d\ — g?2 1/(^1 + c^). Such a system exhibits a 
non-sinusoidal current-phase relation, and a power-law 
decay of the supercurrent with thickness and exchange 
field. If the thicknesses of the both domains are the 
same, the effective field vanishes. Disorder, considered 
both as geometrical fluctuations of the thickness, or ran- 
domly positioned impurities, restores exponential decay 
and sinusoidal phase dependence of the supercurrent. A 
system with two domains of the same width is never in 
the 7r-state. 

To obtain these results, we made a number of sim- 
plifying assumptions. The superconductor- ferromagnet 
interfaces, as well as the boundary between the two fer- 
romagnetic domains, are assumed to be ideal (no scat- 
tering) and sharp. This can be realized in multilayered 
structures, where the ferromagnetic layers can be artifi- 
cially constructed and kept very clean. Another, more 
attractive option, is real ferromagnetic domains. A do- 
main wall has a finite width, typically of order of the 
mean free path, or wider. This induces reflection of elec- 
trons from the domain wall, and additionally generates 
the triplet pairing between electrons and holes. These 
factors need to be taken into account for a quantitative 
comparison between theory and experiment. However, 
we do not expect them to add qualitatively new features 
into the picture we presented. 



2\ 3/2 



J = J0,diff 



x exp 





d I sin I \ j d sin 05. (28) 

hD IV hD ' 



Thus, comparing Eq. I|28l) with Eq. I|27ll we see that a 
long diffusive SFS contact can be a 7r-junction, depending 
on the thickness of the ferromagnet, whereas a similar 
symmetric SFFS contact with anti-parallel configuration 
of the domains is not a 7r-junction. 
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